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ABSTRACT 

Galactic jets are powerful energy sources reheating the intra-cluster medium in galaxy clusters. 
Their crucial role in the cosmic puzzle, motivated by observations, has been established by a great 
number of numerical simulations missing the relativist ic nature of these jets. We present the first 
relativistic simulations of the very long term evolution of realistic galactic jets. Unexpectedly, our 
results show no buoyant bubbles, but large cocoon regions compatible with the observed X-ray cavities. 
The reheating is more efficient and faster than in previous scenarios, and it is produced by the shock 
wave driven by the jet, that survives for several hundreds of Myrs. Therefore, the X-ray cavities in 
clusters produced by powerful relativistic jets would remain confined by weak shocks for extremely 
long periods, whose detection could be an observational challenge. 

Subject headings: Galaxies: active — Galaxies: jets — Hydrodynamics — Shock- waves — Relativistic 
processes — X-rays: galaxies: clusters 



1. INTRODUCTION 

Galaxy clusters are formed by dark matter and gas. 
This last component appears in the form of galaxies and 
a diffuse hot gas filling the space amid them - the in- 
tra cluster medium (ICM). The basic laws of physics 
predict that huge amounts of this gaseous component 
- ICM - should cool due to bremsstrahlung radiation 
and fall onto the central galaxies in the clusters. These 
flows of cold gas would eventually be intimately related 
with crucial processes in the galaxy formation, like for 
instance, the star formation history in galaxies. How- 
ever, these flows - the so called cooling flows - have not 
been observed in most of the clusters, or when observed, 
they are not as important as expected. In order to rec- 
oncile the theoretical results with the observations, sev- 
eral physical mechanisms have been invoked, being the 
most widely accepted, the reheating of th e ICM by the 
Active Galactic Nuclei (A GN) feedback (jFabianI flQOl 
iMcNamara fc NulsenI 120071 ). 

The state-of-the-art picture of the AGN feedback is 
underpinned on the idea that galactic jets can trans- 
port energy from the very center of the galaxy out to 
the cluster scales. These energy injections would in- 
flate bubbles whose evolution would have two well dif- 
ferentiated phases: first, a shock dominated supersonic 
phase, followed by a second subsonic phase when the 
bubbles inflated by th e jets would be buoyant in the ICM 
(jCattaneo et a l."2009^ 

The buoyant bubbles are unstable, due to Rayleigh- 
Taylor and Kelvin-Helmoltz instabilities, when interact- 
ing with the surrounding ICM, leading to a mixture 
of the hot gas locked in the bubble with the environ- 
ment. Besides, an additional mixing is produced at 
the turbulent wake created by the buoyant bubbles ris- 
ing up in the cluster potential well. All this mixing 
produces a net gain of internal energy of the ICM re- 
sulting in an efficient feedback mechani sm able to stop 
or d e lay the cooling flows ( see, e.g., 'Churazov et al.' 
2001[ iQuilis. B ower fc Balogh 2001; Briissen & Kaiser 
20021 : iChurazov et al.l 120021 : IRuszkowski k Begelman I 



2002"; Birzan et al. 
Rovchowdhurv et a 
Sternberg fc Sokerl 



m04; Daha Vecchia et al. 2004; 

^ 12004; Brighenti fc Mathews 2006 ; 

[20081 : iDe Yound I20M The idea 
that the radio lobes reach pressure equilibrium with 
their environment in a relatively short time supports 
the model of buoyant motion of the bubbles inflated by 
the jets (McNamara & Nulsen 2007). However, the ob- 
servations of shocks in several sources like Hercules A 
(jNulsen et al.l [20051 Hvdra A (iSimionescu et alll2009[ ). 
MS0735.6+ 7421 (jMcNamara et al.l [2005D or HCG 62 
(jGitti et al.l l2010l ) may require to reconsider the rele- 
vance of the subsonic buoyant phase. Specially, when 
obser vational evidences (^Kraft et al." 2007:'Croston e t all 
20071) and numerical simulations (Peruch o & Marti 20C 



Bordas. Bosch- Ramon fc Perucholl201l[ ) have shown that 
even very modest energy injections - low-power FRI 
jets - with ages ::^ 10^ yrs still present relatively strong 
shocks. 

The AGN feedback scenario has mainly been stablished 
by a great number of numerical simul ations studying 
the lo ng term evolutio n of jets (e.g., "Churazov et al.| 
20011 : IQuilis. Bower fc Balogh 2001; Briiggen & Kaiser 
2002 |: LRevnolds, Heinz & Begelman 20 02i: lOmma et all 



2OO4I: Omma & Binnev 2004; Zan ni et alJ 12005 



Brighenti & Mathews 2006; Vernaleo & Revnolds 2006; 
Cattaneo fc Tevssier 20()7: Vernaleo fc Revnolds 2007: 



Briiggen e t al.l 12007: iBinnev. Alouani Bibi fc Ominal 
■2007. : .Briiggen et all 120091 : lO^Neih fc JonesI I2010D . In 



these works, the input of jets was modeled injecting 
mass, momentum and energy in a few computational 
cells with a huge size (for jet scales) and with low 
(i.e., non-relativistic) flow velocities and temperatures 
imposed by the Newtonian approach. These inherent 
constrains could have direct implications on the evo- 
lution of the simulated jets since, in order to match 
the typical momentum and energy fluxes, the jets are 
set up with unrealistic opening angles, radii, and flux 
masses when compared with observations. This could 
be the reason leading, in general, to the formation of 
weak shock waves (with Mach numbers Mg < 5) and, 
as a consequence, to an early transition to the subsonic 
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phase in almost all the simulations performed until now. 
In this paper we present, for the first time, the results 
of a set of very long term axisymmetric FRII-like jet sim- 
ulations evolved using a fully relativistic description of 
the fluid dynamics and thermodynamics, also including 
a relativistic equation of state. Our approach allows the 
simulations to reconcile the inferred momentum and en- 
ergy fluxes of observed jets with reasonable values of the 
jet flow velocities, radii, opening angles and mass fluxes. 
The use of the relativistic equation of state accounts con- 
sistently for the relativistic character of electrons in the 
jet and, to a less extent, in the cocoon, and the non- 
relativist ic behaviour of protons. Also our relativistic 
models correctly describe, by construction, the decelera- 
tion of the relativistic flow and the internal and kinetic 
energy fluxes accross the jet terminal shock, that govern 
the eflFective energy flux into the cocoon. 

2. SIMULATIONS 
2.1. Setup 

We simulate the jets as energy injections in a realis- 
tic environment at a distance as close as 1 kpc to the 
nucleus of the source and follow the evolution up to 
several hundreds of kiloparsecs. The jets are injected 
with an initial radius of 100 pc with flow velocities rang- 
ing from Vj = 0.9 c to Vj = 0.99 c, and typical den- 
sity ratio between the jet material and environment of 
Pjl pa = 5 X 10~^. The jets are fed by the injection of 
energy during 50 Myrs (16 Myrs in one model), when 
they are switched off so as to mimic a duty cycle event. 
The total injected energies range from 3 x 10^^ to 10^^ 
erg, depending on the jet model (see solid-black lines in 
Fig. [3]). Our best numerical resolution is 50 x 50 par- 
sees. A complete list of all relevant parameters of the 
simulations is described in Table 1. 

All the performed simulations are 2D axisymmetric. 
The jets are injected in a computational domain filled 
with an ambient in hydrostatic equilibrium, which is 
formed by primordial gas with a King-like density pro- 
file considering the elliptical galaxy - origin of the jet 
- and the galaxy cluster. The density profile param- 
eters have be en fixed by fitting the X-ray data of the 
source 3C 31 fHardc astle et al.ll2002[ ). The ambient den- 
sity at the injection point is 0.1 particles per cubic cen- 
timetre . The dark matter hal o follows a NEW density 
profile (|Navarro. Frenk fc Whi te 1997). All these param- 
eters represent a moderate size galaxy cluster with mass 
10^^ Mq and ~ 1 Mpc virial radius. 

The numerical grid is structured as follows: in the ra- 
dial direction, a grid with the finest resolution extends 
up to 50 kpc (Model Jl) or 100 kpc (Models J2, J3, 
J4). An extended grid with decreasing resolution was 
added up to 1 Mpc. The time-step during the first part 
of the simulations, when the jet is still active, was 50 to 
100 years. The boundary conditions in the simulations 
are reflection at the jet base, to mimic the presence of a 
counter-jet, reflection at the jet axis and outflow at the 
end of the grid in the axial and radial directions. 

The simulations presented in this paper use the finite- 
volume code Ratpenat^ which solves the equations of rela- 
tivistic hydrodynamics in conservation form using high- 
resolution-shock-capturing methods. Ratpenat is a hy- 
brid parallel code - MPI + OpenMP - extensively and 



intensively tested (e.g., iPerucho et al.l 120101) . The code 
includes the Synge equation of state ()Svngeiri957[ ) with 
two populations of particles, namely, leptons (electrons 
and positrons) and baryons (protons). In these simula- 
tions, cooling has not been taken into account, as the typ- 
ical cooling times in the environment a re long compared 
to the simulation times (see Figure 10 in lHardcastle et aTl 
^002,). 

2.2. Results 

The dynamics of the system is dominated by the jet 
active phase, where jet head velocities range from 0.01 c 
to 0.06 c, consistently with different estimates of the 
advance velocities of ac tive radio sources (0.02 — 0.2 c) 
(jCarilli fc Bart hell Il996f ). implying Mach numbers be- 
tween 10 (J3) and 30(J2). The evolution of the super- 
sonic jet generates a characteristic morphology: i) a bow- 
shock that acts on the ambient medium, ii) a terminal 
or reverse shock at the head of the jet where the flow is 
decelerated and iii) the cocoon inflated by the shocked 
jet particles and polluted with shocked ICM stirred via 
instabilities arising at the contact discontinuity between 
both media, typically hotter and underdense compared 
with the ambient. After the switching off, due to the 
short time scales needed by the relativistic fluid to reach 
the jet terminal shock, the jet head velocities quickly 
drop to values ~ 10~^ c making the bow-shock to ap- 
proach sphericity very fast. During this phase, the Mach 
numbers of the bow-shock fall from ~ 10 to values be- 
tween 1 and 2. 

Figure 1 shows four snapshots of model J2 at represen- 
tative phases of its evolution (see Sec. 3.1). The jet is 
seen in the second panel (also in the first one although 
less clearly due to the small size of the system) as a blue 
or green line on the axis in the density and temperature 
frames, respectively. The terminal shock at the head of 
the jet is also seen in this panel, specially in the temper- 
ature distribution, as the dark red (saturated) region at 
the end of the jet. After the jet switch-off, the channel 
opened by the jet is still seen on the axis in the third and 
fourth frames although the jet terminal shock has already 
disappeared. In the four panels, the cocoon is the bluish 
turbulent region in the density frames (reddish region 
in the temperature ones), whereas the shocked ambient 
medium fills the region between the cocoon and the bow 
shock (yellow/orange/red region in the density frames; 
light blue region in the temperature frames). The den- 
sity is low (i.e., smaller than the density in the origi- 
nal unperturbed medium) in the cocoon and high in the 
shocked ambient region. 

The key features that differentiate the simulated jets 
are the injection power and the duration of the active 
phase. The location of the bow shock at a given time 
is largely dependent on these parameters and hence the 
distance at which the energy is deposited. The morphol- 
ogy of the cocoon also changes for the different models, 
with the most powerful jet (J2) creating a large cavity 
and the less powerful one (J3) having its particles dis- 
tributed closer to the source. The gross properties of 
models Jl and J4, with the same power and different 
composition, are very similar. 

In all the cases studied, the late stages of the simula- 
tions show a cocoon composed by jet and ICM shocked 
particles completely mixed, forming a low density region 
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TABLE 1 

Parameters of the simulated jets 



Model 


Velocity [c] 


Density [g/cm^] 


Xe 


Lfc [erg/s] 


max. resol. [pc/cell] 


toff [Myrs] 


Jl 


0.9 


8.3 X 10"^^ 


1.0 


10^^ 


50 


50 


J2 


0.984 


8.3 X 10-"^ 


1.0 


^q46 


100 


16 


J3 


0.9 


8.3 X 10-^° 


1.0 


^q44 


100 


50 


J4 


0.9 


8.3 X 10-2^ 


0.5 


10^^ 


100 


50 



From left to right the columns give the model, the injection velocity, the injection density, the leptonic number, the jet power, the maximum 
resolution, and the switch-off time. 



surrounded by the denser shocked ICM, with tempera- 
tures within the cocoon not larger than one order of mag- 
nitude those in the shocked medium. All these features 
postulate the cocoon as an excellent candidate precur- 
sor of the X-ray observed cavities. Surprisingly, in all 
our simulations, covering two orders of magnitude in jet 
power and a factor of 30 in total injected energy, the pres- 
sure jumps between the shocked and the unperturbed 
ICM persist along the whole simulation, and no buoyant 
stage is reached. Fig. 2 shows the X-ray emission ob- 
tained from one snapshot in simulation J2 compared to 
a real observation. 

The minimum cooling times have been calculated at 
different times during the simulation to check the consis- 
tency of the results, taking into account that cooling has 
not been included. After the medium is perturbed by 
the jet, minimum cooling times are longer than the sim- 
ulated ones by more than a factor 10 during the active 
phase, and by more than 100 or 1000 times the simulated 
time after the jet switch-off. Significant changes in the 
results presented in this work are thus not expected if 
the cooling terms are included. 

3. DISCUSSION 
3.1. Cavity evolution 

The long-term evolution of the cocoons in our numer- 
ical simulations, within the shock dominated supersonic 
phase, can be interpreted as undertaking three stages: i) 
a short one-dimensional phase, in which the cocoon ex- 
pansion is governed by the one-dimensional evolution of 
the jet (up to t ^ 1.8 Myrs in model J2); ii) a genuinely 
multidimensional phase in which the cocoon expansion 
is driven by a decelerating jet as a result of the multi- 
dimensional effects affecting the jet propagation, and iii) 
a Sedov phase (starting at t = 16 Myrs in model J2), 
in which the cocoons expand passively once the energy 
injection of the jet has ceased. In this last phase, the ex- 
pansion on a density decreasing atmosphere is expected 
to produce a faster expansion of the cocoon (and conse- 
quently a faster pressure decrease) than in the pure Sedov 
case (constant ambient density). The first two phases are 
typical of the propagation of supersonic jets and are de- 
scribed in, e.g., Scheck et al. (2002). The third phase is 
new and follows from the jet switch-off. 

Figure 1 shows four snapshots of model J2 at rep- 
resentative stages of its evolution: one-dimensional 
(first panel), multidimensional (second panel) and Sedov 
phases (third and fourth panels). In a homogeneous am- 
bient, during the so-called one-dimensional phase, the 
jet propagates at its estimated one-dimensional speed 
(iMarti et al.iil997i ). The propagation efficiency of the 
jet during this phase is maximum and the cocoon in- 



flates at its smallest rate. In the present simulations, the 
evolution in a density decreasing atmosphere accelerates 
the jet propagation speed beyond the one-dimensional 
estimate and makes the cocoon expansion faster. The 
onset of the multidimensional phase is triggered by mul- 
tidimensional, dynamical processes taking place close to 
the jet's head. During this phase, the generation of large 
vortices at the jet/shocked ambient interface (a pair of 
such vortices are seen at the jet's head in the density 
map of the second panel) decelerates the jet increasing 
the flux of jet material into the cocoon. As a result, the 
cocoon expands at a fast rate (helped also by the de- 
creasing ambient density). Third and fourth panels in 
Fig. 1 are representative of the Sedov phase, once the 
jet has been switched-off. Consistently, the jet losses 
its hot-spot (the impact point of the jet on the ambient 
medium), still seen in the second panel, and the chan- 
nel opened by the jet starts to be refilled. During this 
phase, the cocoon continues its expansion (at a smaller 
rate than in the previous phase since the injection of en- 
ergy has ceased) and tends to sphericity (since there is 
no extra momentum transfer in the axial direction any 
longer). 

The temperature is quite homogeneous within the co- 
coon, and almost constant with time as far as the there 
is a continuous and constant matter and energy injection 
through the jet. Once the jet is switched-off, the cocoon 
starts to cool down (see the evolution of the cocoon tem- 
perature between the third and fourth panels). 

At late stages, the velocity field in the cocoon is 
dominated by the overall cocoon expansion speed (< 
3 X 10~^ c) and the turbulent motions on smaller scales, 
with local values of the order of 0.04 c. 

Begelman-Cioffi's model (|Begelman fc Cio"ffil Il989f ) 
(BC) describes the expansion against a uniform ambi- 
ent medium of the overpressured cocoons raised by the 
continuous injection of energy from a supersonic jet. In 
this model the axial expansion of the cocoon (i.e. along 
the jet) proceeds at the (constant) advance speed deter- 
mined by the jet, whereas the sideways growth follows 
from the assumption of the evolution being mediated by 
a strong shock. 

The long-term evolution of the cocoons in our 
numerical simulations can be consistently described 
withi n the so-called extended B egelman-Cioffi's model 
(Sche ck et alllM)! [Perucho fc Marti 2007) (eBC) that 
allows for a power-law dependence of the jet advance 
speed with time and a non-uniform ambient medium. In 
addition, the model can also describe the passive (super- 
sonic) expansion of the cocoon once the jet has ceased 
its activity (Sedov phase). 

In the eBC model, if the advance speed of the bow 
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Fig. 1. — Maps of logarithms of rest-mass density (left) and temperature (right) for simulation J2 at times 1.1, 13.5, 34 amd 180 Myrs. 
The figures shows a mirrored image around the symmetry axis and plane in the simulation. The color-scale (in the online version only) of 
the temperature plot has been cut at 10-'^'-' K for the sake of clarity. Although the maxima in the first and second snapshots are 5 x 10-'^'-' K 
and 4.5 X 10-*^ -"^ K, this values are only reached in a very small region at the head of the jet (hotspots), which accordingly saturate the 
color-scale used. 



shock along the axial direction, Vc 
sity, Pa, follow the power laws Vc 



and the ambient den- 



(X t"^, Pa CX 



,/5 



then, 



the cavity's transversal dimension. Re, and pressure, Pc, 
follow 



Re (xt 4+? 



2(c.-2)-c.(4+/3) 
Pc (Xt 4+/3 



(1) 



These time dependencies are valid as far as there is a 
constant injection of energy in the cocoon. Once the jet 
is switched off (Sedov phase), the evolution follows 



Re (Xt^+f3 , Pc (X t~ 



2(a-l)-(l + a)(4+/3) 



(2) 



(the pure Sedov expansion phase is recovered when a = 
-3/5 and /3 = 0, for which R^ ex t^/^ and P^ ex t'^/^). 

Table 2 displays the parameters a (derived from the 
simulations) and /3 for the three phases of each simula- 
tion, and the corresponding exponents for the time de- 



pendence of Re and Pc both from the simulations and 
from the eBC model. According to the data displayed in 
this table, we can conclude that the eBC model describes 
consistently the long-term evolution of the cocoons along 
phases i) and ii). Concerning the Sedov phase, let us 
note first that our model produces a faster expansion 
(time exponent for Re around 3/5 instead of 2/5) and 
a faster pressure decrease (time exponent for Pc around 
—7/5 instead of —6/5), as expected. However it must be 
noted that the exponents derived from the simulations 
reinforce this tendency. The discrepancies between the 
expected and the obtained time dependencies in this last 
phase have a clear dependence on jet power and could 
be a signature of the buoyant force acting on the (still 
supersonic) cavities. 
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TABLE 2 

Parameters for the three stages of the evolution within each of the simulations. 





a 


ID 


phase 

Pc 


Re 


a 


2D 


phase 

Pc 


Re 


a 


Sedov 


phase 

Pc 


Re 


Jl Sim 
Model 


0.07 


-1.55 


-1.58 
-1.65 


0.75 
0.79 


-0.23 


-0.52 


-1.09 
-1.05 


0.66 
0.64 


-0.74 


-1.02 


-1.70 
-1.43 


0.90 
0.58 


J2 Sim 
Model 


0.27 


-1.55 


-1.67 
-1.68 


0.67 
0.71 


-0.57 


-0.52 


-0.95 
-0.91 


0.81 
0.74 


-0.83 


-1.02 


-1.67 
-1.40 


0.72 
0.61 


J3 Sim 
Model 


0.13 


-1.55 


-1.55 
-1.66 


0.67 
0.76 


-0.35 


-0.52 


-1.08 
-1.00 


0.74 
0.68 


-0.60 


-1.02 


-2.16 
-1.47 


1.00 
0.54 




Fig. 2. — Qualitati ve comparison of a n X-ray map of the clus- 
ter MS 0735.6+7421 ([McNamara et al.i r2005) (left, credit: X-ray 
image: NASA / CXC / Ohio U. / B.McNamara et al.; illustra- 
tion: NASA / CXC / M.Weiss) and a synthetic X-ray luminosity 
map extracted from simulation J2 (right). Although this compari- 
son must be taken with caution, the main features of both images 
seem to match remarkably well. 



The parameters ol and (3 are derived from the simulations (first two columns for each phase). The time dependence of Re and Pc (third 
and fourth columns, respectively) is shown as obtained from the simulation and from the eBC model. 

instabilities arising at the contact discontinuity between 
the shocked ambient and the shocked jet fluid. There- 
fore, the bulk of the heating comes from the action of the 
shock, being the mixing with the jet material negligible 
from the energetic point of view. 

Figure [H shows the energy per particle inside the 
shocked region versus time for all four simulations. The 
energy per particle, including the jet and ambient com- 
ponents, inside the bow-shock is always over 1 keV for 
Jl, J2 and J4, whereas it falls below this value after 
c^ 10^ yrs in J3 case. Any effective reheating mechanism 
able to stop the cooling flows would require energies per 
particle r^ 1 keV (McNamara fc Nulsen. 2007). Figure [5] 
shows the average internal energy per unit volume ver- 
sus distance to the source at different times, for the case 
of J2. The plot reveals where the energy is deposited. 
The values at small radii at the smaller times (mainly 
1 and 15 Myr) are influenced by the elongated shape of 
the bow-shock. When the bow-shock becomes closer to 
sphericity at later times, the energy density has its max- 
imum in a wide shell behind the shock- front. These facts 
confirm that the bow-shock dominates the energy deposi- 
tion during the whole evolution. The heating of the cool- 
ing region (::^ 20 kpc, Hardcastle et al. 2002) is very fast 
(< 1 Myr). It is also remarkable that, after 50 Myrs, the 
energy density falls below the original one. In the case 
of J2, the injection time was fixed to be 16 Myrs, imply- 
ing that the simulation has been run for c:± 12 times the 
active time. Taking into account that the heating mech- 
anism is very fast and acts mainly at the bow-shock, as 
the shock propagates outwards, the central regions cool 
down due to expansion. This fast cooling is the rea- 
son why sl ower heating mechanism s are favoured in the 
literature rtOmma fc Binnevll2QQ4f ). Nevertheless, there 
are two aspects to be taken into account: 1) the gas 
in the inner region is still expanding, and 2) consider- 
ing that this activity period has only implied the injec- 
tion of :^ 2 X 10^ M0, which does not represent a large 
amount of material out of the total amount of matter 
and energy expected to be present in the surroundings 
of AGN. New periods of activity would be expected to 
occur more often than 16 Myrs every 180 Myrs (see, e.g., 
iFabian eFal. 2006) . We note that even short periods of 
activity (< 1 Myr) with the same power, are enough to 
heat and empty the ambient inside the cooling region. 

The moderate size galaxy cluster considered in our 
simulations, with mass 10^^ Mq and 1 Mpc virial radius, 
does favour the persistence of the bow shock generated in 
the expansion of the cavities. However, a ten times more 



3.2. Heating of the ambient medium 

The efficiency of the heating triggered by the bow 
shock is very high in two aspects, namely the amount of 
ambient medium heated per unit time and the amount of 
energy transferred to the ambient medium in the heat- 
ing process. With respect to the first point, let us note 
that the heating front propagates through the ambient 
medium at the shock speed (a substantial fraction of 
the light speed during the first stages of evolution). As 
an example, the shock in model J2 has processed about 
9.4 X 10"^ kpc^ by the end of the one-dimensional phase 
{t ^ 1.8 Myrs), 1.1 x 10^ kpc^ at the onset of the Sedov 
phase {t ^ 16 Myrs), and 0.37 Mpc^ at the end of the 
simulation {t ~ 160 Myrs). 

The transfer of energy to the ambient medium is also 
very efficient. Indeed, the energetic balance of our simu- 
lations (Fig. 3) shows that between 95% and 97% of the 
injected energy by the jet is instantaneously transferred 
to the ambient medium through shock-heating, mixing, 
and acceleration. The small residual is invested on gain- 
ing potential energy or kept by the jet particles. By 
the end of our simulations, between 10^^ to 10^^ Mq 
of ICM gas have been heated up by the shock, depend- 
ing on the total energy injected during the active phase, 
in ag reement with constrains ini posed by observational 
data ("McNamara fc N ulsen|[2QQ7h . Only a small fraction 
of the reheated ICM gas - from 0.1 to 1% - is mixed by 
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Fig. 3. — Logarithm of energy versus time for the four simulations. The top panels show the results for models Jl and J2, respectively, 
whereas the bottom panels show those for J3 and J4. The red-dotted lines (color only in the online version) and the red-dashed lines 
represent, respectively, the increase of internal and kinetic energy in the processed ambien medium. The blue lines, dotted and dashed, 
display the internal and kinetic energy of the jet material. The increase of potential energy is shown by the black dashed-dotted line. The 
upper thick black line is the total injected energy. These plots reveal that the jet material barely keeps a small fraction of its injected 
energy (less than 1%), which is mainly transferred to the ambient. After the switching off, the jet material keeps on transferring energy 
to the ambient via mixing, being this the explanation to the energy drop after this time. A tracer, evolved in the code as an additional 
conserved variable in the set of equations, allow us to accurately discriminate between jet and ambient material. 
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Fig. 4. — Energy per particle inside the bow-shock versus time 
for the four simulations. Jl is represented by red color, J2 by black, 
J3 by yellow, and J4 by blue (color only in the online version). 



massive cluster with an average density (and pressure) 



10 times larger will reduce the volume processed by the 
shock in the same factor or, equivalently, in roughly a fac- 
tor of 2 per spatial dimension. Moreover large changes 
in the processed mass are not expected because of the 
denser ambient medium. Finally, the same amount of 
injected energy acting on a similar amount of particles 
would lead to similar values of energy per particle. 

3.3. Jet /cavity energy balance 

The energy injected by the jet is invested primarily 
in inflating the cavity against the ambient medium and 
heating the ambient gas remaining within it. Other con- 
tributions, like, e.g., the kinetic energy transferred to the 
ambient, the change in the potential gravitational energy 
of the pushed-aside material, and the energy kept by the 
jet particles, are neglected in the jet/cavity energy bal- 
ance estimates in real sources. The energy invested in in- 
flating the cavity is the p F expansion work exerted on the 
ambient. It is usually estimated from the volume of the 
cavity and the average pressure surro unding it. In the 
case of the cluster MS 0735.6+7421 ( McNamara et al.l 
l2005f ). where the cavities are roughly 200 kpc in diame- 



ICM reheating by relativistic jets. 
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Fig. 5. — Average internal energy density versus radius from the 
nucleus for simulation J2, at differet times. The blue-solid line 
is calculated at t = 1 Myr, the green-dotted line at t = 15Myr, 
the blue-dashed line at f = 50 Myr, and the red-solid line at t = 
180 Myr (color only in the online version). 



ter and the surrounding pressure, 6 x 10~^^ ergs/cm^, the 
work required to inflate each cavity against this pressure 
is pV ~ 10^^ ergs. The internal energy in the cavity, i.e., 
the second term in the jet/cavity energy balance, is usu- 
ally estimated as being few times the pV term (3 times 
in the case of the MS 0735.6+7421 cluster). 

Our simulations allow for a check of the hypothesis be- 
low the jet/cavity energy balance estimates summarized 
in the previous paragraph. First of all the results shown 
in Fig. 3 confirm that the kinetic energy in the ambi- 
ent gas, the change in the potential gravitational energy, 
and the energy kept by the jet particles are negligible, 
contributing with less than 10% to the global energetics. 
Now, the expansion work of the cavity can be approx- 
imately accounted for through estimates of the cocoon 
volume and the average pressure of the original ambient 
medium swept by the jet. For the model J2 {t = 180 
Myrs), the volume of the cocoon is ~ 3.4 x 10^ kpc^ 
(i.e., about 400 kpc of effective diameter), whereas the 
averaged pressure is ~ 3 x 10~^^ erg/cm^ (200 times the 
pressure used in the MS 0735.6+7241 calculation). The 
resulting work is pV ~ 3 x 10^^ ergs. Finally, the av- 
erage internal energy density in the cavity is estimated 
to be 10~^^ erg/cm^. With these numbers, the mini- 
mum total energy needed to inflate one of the cavities in 
simulation J2 results to be 1.3 x 10^^ ergs, in good agree- 
ment with the total energy injected by the jet (5 x 10^^ 
ergs), reflecting the consistency of our analysis. Finally, 
as a by-product, we note that our result confirms the 
factor of 3 between the internal energy stored in the cav- 
ity and the expansion work assumed in the case of the 
MS 0735.6+7421 cluster, but now under very different 
conditions: in a two orders of magnitude lighter cluster 
and for a total injected energy 10 times smaller. 

3.4. Possible three-dimensional effects 

We must make a final comment on the 2D nature of the 
simulations, for it could be seen as an important limita- 
tion in our conclusions. In the long term, the evolution of 



the cavity is driven by its internal pressure, which i) de- 
pends basically of the total energy injected in the cocoon, 
and ii) tends to isotropize the cocoon on a characteristic 
time scale of its internal sound speed (of the same order 
or larger than the overall shock advance speed). Both 
aspects tend to reduce the importance of the early three- 
dimensional effects in the long term evolution of cavities 
making more reliable the conclusions derived from our 
present axisymmetric simulations. Moreover, the ambi- 
guity in the initial condition parameter space for the jet 
and ambient could lead to results with larger dispersion 
than those arising from possible 3D effects. 



4. CONCLUSIONS 

We have presented the first simulations on the impact 
of relativistic AGN jets in the heating of the intracluster 
medium. Our simulations cover the longest spatial and 
temporal scales considered up to now and besides this, 
due to their relativistic character, also for the first time, 
allow for a consistent description of the jets (that have 
consistent values of the jet flow velocities, radii, opening 
angles and mass, momentum and energy fluxes) within 
the cluster heating scenario. 

Our simulations show that heating by AGN jets is 
mainly driven by shock heating, resulting in a very fast 
and efficient process (more than 95% of the energy in- 
jected through the jets is transferred to the ambient). 
As a by-product, our simulations also show that although 
buoyant cavities could be the very final stage of radio- 
galactic relics, they are not the main actors in the ICM 
reheating. 

As a result of our simulations, we suggest the idea 
that most of the observed X-ray cavities are confined 
by shock waves, very weak though. Although it can not 
be concluded that the presence of shock confined cav- 
ities implies the relativistic nature of jets, such shocks 
in our results last for much longer periods than in pre- 
vious - Newtonian - works. The conffimation of the 
existence of such weak shock waves could be an observa- 
tional challenge that would have crucial implications in 
our understanding of the galaxy formation and evolution 
paradigm. 



Simulations were carried out using resources from the 
Spanish Supercomputing Network (RES). This Research 
Project has been supported by Spanish Ministerio de 
Ciencia e Innovacion (grants AYA2010-21322-C03-01, 
AYA2010-21097-C03-01 and CONSOLIDER2007-00050) 
and by the Generalitat Valenciana (grant PROMETEO- 
2009-103). MP acknowledges support from MICINN 
through a "Juan de la Cierva" contract. We thank B.R. 
McNamara for providing an X-ray image of the cluster 
MS 0735.6+7421. The authors thank the anonymous ref- 
eree for constructive criticism and comments that have 
improved the manuscript. 



Perucho, Quilis and Marti 



REFERENCES 



Begelman, M.C., Cioffi, D.F., 1989, ApJ, 345, L21 

Binney, J., Alouani Bibi, P., Omma, H., 2007, MNRAS, 377, 142 

Birzan, L., Rafferty, D.A., McNamara, B.R., Wise, M.W., Nulsen 

P.E.J. , 2004, ApJ, 607, 800 
Bordas, P., Bosch-Ramon, V., Perucho, M., 2011, MNRAS, 412, 

1229 
Brighenti, P., Mathews, W.G., 2006, ApJ, 643, 120 
Briiggen, M., Kaiser, C.R., 2002, Nature, 418, 301 
Briiggen, M., Heinz, S., Roediger, E., Ruszkowski, M., 

Simionescu, A., 2007, MNRAS, 380, L67 
Briiggen, M., Scannapieco, E., Heinz, S., 2009, MNRAS, 395, 2210 
Cariin, C.L., Barthel, P.D., 1996, A&A Rev., 7, 1 
Cattaneo, A., Teyssier, R., 2007, MNRAS, 376, 1547 
Cattaneo, A., Paber, S.M., Binney, J., et a/., 2009, Nature, 460, 

213 
Churazov, E., Briiggen, M., Kaiser, C.R., Bohringer, H., Porman, 

W., 2001, ApJ, 554, 261 
Churazov, E., Sunyaev, R., Porman, W., Bohringer, H., 2002, 

MNRAS, 332, 729 
Croston, J.H., Kraft R.P., Hardcastle M.J., 2007, ApJ, 660, 191 
Dalla Vecchia, C, Bower, R.G., Theuns, T., Balogh, M.L., 

Mazzotta, P., Prenk, C.S., 2004, MNRAS, 355, 995 
De Young, D.S., 2010, ApJ, 710, 743 
Pabian, A.C., 1994, ARA&A, 32, 277 
Pabian, A.C., Sanders, J.S., Taylor, G.B., Allen, S.W., Crawford, 

C.S., Johnstone, R.M., Iwasawa, K., 2006, MNRAS, 366, 417 
Gitti, M., O'Sullivan, E., Giacintucci, S., David, L.P., Vrtilek, 

Jan, Raychaudhury, S., Nulsen, P.E.J. , 2010, ApJ, 714, 758 
Hardcastle M. J., Worrall D. M., Birkinshaw, M., Laing, R. A., 

Bridle, A. H. 2002, MNRAS, 334, 182 
Kraft, R.P., Nulsen, P.E.J., Birkinshaw, M., 2007, ApJ, 665, 1129 
Marti, J.M., Miiller, E., Pont, J. A., Ibaiiez, J.M., Marquina, A. 

1997, ApJ, 479, 151 



McNamara, B.R., Nulsen, P.E.J. 2007, ARA&A, 45, 117 
McNamara, B.R., Nulsen, P.E.J. , Wise, M.W. et al., 2005, 

Nature, 433, 45 
Navarro, J. P., Prenk, C.S., White, S.D.M., 1997, ApJ, 490, 493 
Nulsen, P.E.J. , Hambrick, D.C., McNamara, B.R. et al., 2005, 

ApJ, 625, L9 
Omma, H., Binney, J., Bryan, G., Slyz, A., 2004, MNRAS, 348, 

1105 
Omma, H., Binney, J., 2004, MNRAS, 350, L13 
O'Neill, S.M., Jones, T.W., 2010, ApJ, 710, 180 
Perucho, M., Marti, J.M., 2007, MNRAS, 382, 526 
Perucho, M., Marti, J.M., Cela, J.M., Hanasz, M., de la Cruz, R., 

Rubio, P., 2010, A&A, 519, A41 
Quilis, v.. Bower, R.G., Balogh, M.L., 2001, MNRAS, 328, 1091 
Reynolds, C.S., Heinz, S., Begelman, M.C., 2002, MNRAS, 332, 

271 
Roychowdhury, S., Ruszkowski, M., Nath, B.B., Begelman, M.C., 

2004, ApJ, 615, 681 
Ruszkowski, M., Begelman, M.C., 2002, ApJ, 581, 223 
Scheck L., Aloy, M.A., Marti, J.M., Gomez, J.L., Miiller, E., 

2002, MNRAS, 331, 615 
Simionescu, A., Roediger, E., Nulsen, P.E.J, et al. 2009, A&A, 

495, 721 
Sternberg, A., Soker, N., 2008, MNRAS, 389, L13 
Synge, J.L., 1957, The Relativistic Gas. North-Holland publishing 

Company, Amsterdam 
Vernaleo, J.C, Reynolds, C.S., 2006, ApJ, 645, 83 
Vernaleo, J.C, Reynolds, C.S., 2007, ApJ, 671, 171 
Zanni, C, Murante, G., Bodo, G., Massaglia, S., Rossi, P., 

Perrari, A., 2005, A&A, 429, 399 



